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Abstract 

The angular distribution and differential branching fraction of the decay 
B° — > K*°/j, + /j,~ are studied using a data sample, collected by the LHCb experi- 
ment in pp collisions at y's = 7 TeV, corresponding to an integrated luminosity of 
l.Ofb -1 . Several angular observables are measured in bins of the dimuon invariant 
mass squared, q 2 . A first measurement of the zero-crossing point of the forward- 
backward asymmetry of the dimuon system is also presented. The zero-crossing 
point is measured to be q$ = 4.9 ± 0.9 GeV 2 /c 4 , where the uncertainty is the sum of 
statistical and systematic uncertainties. The results are consistent with the Standard 
Model predictions. 
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1 Introduction 



The B° —7- K*°[i + fi~ decayQ where K*° — > K + ir~ , is a b — > s flavour changing neutral 
current process that is mediated by electroweak box and penguin type diagrams in the 
Standard Model (SM). The angular distribution of the K + n~ fi + fi~ system offers particular 
sensitivity to contributions from new particles in extensions to the SM. The differential 
branching fraction of the decay also provides information on the contribution from those 
new particles but typically suffers from larger theoretical uncertainties due to hadronic 
form factors. 

The angular distribution of the decay can be described by three angles (9g, 9k and 
0) and by the invariant mass squared of the dimuon system (q 2 ). The B° — > K* Q ^ + ^i~ 
decay is self-tagging through the charge of the kaon and so there is some freedom in the 
choice of the angular basis that is used to describe the decay. In this paper, the angle 
9g is defined as the angle between the direction of the fi + in the dimuon rest frame 
and the direction of the dimuon in the B° (B°) rest frame. The angle 9k is defined as 
the angle between the direction of the kaon in the K*° (K*°) rest frame and the direction 
of the K*° (K*°) in the B° (B°) rest frame. The angle is the angle between the plane 
containing the /i + and /i~ and the plane containing the kaon and pion from the K*°. A 
detailed description of the angular basis is given in Appendix [Aj In this basis, the angular 
definition for the B° decay is a CP transformation of that for the B° decay. 

Using the notation of Ref. [l], and assuming equal numbers of B° and B° decays, the 
differential decay rate corresponds to 



d 4 r 



dr / dg 2 dg 2 d cos 9g d cos 9k 



32tt 



3{ sin 2 9 K + S$ cos 2 9 K + 

S% sin 2 K cos 29 g + cos 2 9 K cos 29 1 + 

S3 sin 2 9k sin 2 9i cos 20 + S4 sin 29k sin 29 £ cos <p + 

S5 sin 29k sin 9g cos cf> + Sq sin 2 9k cos 9g + 

Sj sin 29k sin 9e sin + Ss sin 29k sin 26 e sin + 

S9 sin 2 9k sin 2 6g sin 20 , 

(1) 



where the S{ terms are CP averages between B° and B° of bilinear combinations of K*° 
decay amplitudes that vary with q 2 . The terms 5*7, Sg and Sg are suppressed by the small 
size of the strong phase difference between the amplitudes involved and are expected to 
be close to zero across the full q 2 range not only in the SM but also in most extensions. 
To reveal the effect of new particles, it is better to look instead at the corresponding CP 
asymmetries A 7 , A 8 and Ag, between B° and B°, which are not suppressed by the size of 



1 Charge conjugation is implied throughout this paper unless stated otherwise. 
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29 the strong phases involved [21. This is achieved by flipping the sign of <fi for the B decay 

30 in Eq. [TJ replacing S 7 , S 8 and S 9 by A 7 , A 8 and A 9 . 

31 If q 2 ^> 4m 2 (where m M is the muon mass), the terms S^, S[, S£ and S| in Eq. [ljare 

32 related to the fraction of longitudinal polarisation of the K*° meson, (S^ = — S£ = Fl 

33 and fSJ = 4Sf = 1 — Fl). The angular term sin 2 9 k cos 9i is responsible for generating 

34 the forward-backward asymmetry of the dimuon system, A FB |3J (A FB = f S 6 ). The term 

35 S3 is proportional to the asymmetry between the two transverse K*° amplitudes, A\ jij. 

36 Due to the limited number of signal candidates that are available in the present analysis, 

37 the angle <fi is also transformed such that 

if ^ <0 (2) 
otherwise 

38 to cancel terms in Eq. [T] that have either a sin or a cos dependence and provide a 

39 simplified angular expression. The reduced angular expression then contains only F^, A f -q, 
S3 and Ag and is given in full in Sec. [7j 

41 In the large-recoil (low q 2 ) limit, the observables A\ and A^ e , where S3 = \ (1 — Fj) A\ 

42 and A F b = § (1 — Fl) ^t 2 ' can be defined and have reduced form-factor uncertainties |4||5]- 

43 The observables A\ and A^ e depend only on the transverse decay amplitudes. 

44 This paper presents a measurement of the differential branching fraction (dB/dq 2 ), 

45 A F b, -Fl, S3 and Ag of the B° — > K*°fi + \T decay in six bins of q 2 . Measurements of 

46 the transverse asymmetries A\ and ^4^ e are a ^ so presented. The analysis is based on a 

47 dataset, corresponding to l.Ofb -1 of integrated luminosity, collected by the LHCb detector 

48 in y/s = 7TeV pp collisions in 2011. Section [2] describes the experimental setup used in 

49 the analyses. Section [3] describes the event selection. Section [4] discusses potential sources 

50 of peaking background. Section [5] describes the treatment of the detector acceptance in 

51 the analysis. Section [6] discusses the measurement of dB/dq 2 . The angular analysis of 

52 the decay, in terms of cos6^, cos 9k and 0, is described in Sec. [7J Finally, in the SM, A F b 

53 varies as a function of q 2 and is known to change sign at a predictable value of q 2 . A first 

54 measurement of this zero-crossing point is presented in Sec. [8] 

55 The one- dimensional cos 9g, cos 9k and distributions have previously been studied at 

56 LHCb |6|, BaBar [7], Belle [8] and CDF [9] with smaller data samples. 



2 The LHCb detector 



58 The LHCb detector 10 is a single-arm forward spectrometer, covering the pseudorapidity 

59 range 2 < r\ < 5, that is designed to study b and c hadron decays. A dipole magnet with a 

60 bending power of 4 Tm and a large area tracking detector provide momentum resolution 

61 ranging from 0.4% for tracks with a momentum of 5GeV/c to 0.6% for a momentum 

62 of lOOGeV/c. A silicon microstrip detector, located around the pp interaction region, 

63 provides excellent separation of B meson decay vertices from the primary pp interaction 

64 and impact parameter resolution of 20 (Jin for tracks with high transverse momentum (pr)- 



65 Two ring-imaging Cherenkov (RICH) detectors 11 provide kaon-pion separation in the 
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67 
68 
69 
70 
71 
72 
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momentum range 2 — lOOGeV/c. Muons are identified based on hits created in a system 
of multiwire proportional chambers interleaved with layers of iron. The LHCb trigger |12 
comprises a hardware trigger and a two-stage software trigger that performs a full event 
reconstruction. 

Samples of simulated events are used to estimate the contribution from specific sources 
of exclusive backgrounds and the efficiency to trigger, reconstruct and select the B° — > 



K*°n + fi~ signal. The simulated pp interactions are generated using Pythia 6.4 13 with 
a specific LHCb configuration 14 . Decays of hadronic particles are then described by 

Finally, the 
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EvtGen 15 in which final state radiation is generated using Photos 
Geant4 toolkit |17] is used to simulate the detector response to the p 
by Pythia/EvtGen, as described in Ref. 18 . The simulated samples are corrected 
for known differences between data and simulation in the B° momentum spectrum, the 
detector impact parameter resolution, particle identification 11 and tracking system 
performance using control samples from the data. 



so 3 Selection of signal candidates 

si The B° — > K*°fi + fi~ candidates are selected from events that have been triggered by a 

82 muon with p^ > 1.5GeV/c, in the hardware trigger. In the first stage of the software 

83 trigger, candidates are selected if there is a reconstructed track in the event with high 

84 impact parameter (> 125 |xm) with respect to one of the primary pp interactions and 

85 p-r > 1.5 GeV/c. In the second stage of the software trigger, candidates are triggered on 
se the kinematic properties of the partially or fully reconstructed B° candidate [l2| . 

s? Signal candidates are then required to pass a set of loose (pre-) selection requirements. 

88 Candidates are selected for further analysis if: the B° decay vertex is separated from 

89 the primary pp interaction; the B° candidate impact parameter is small, and the impact 

90 parameters of the charged kaon, pion and muons are large, with respect to the primary pp 

91 interaction; and the angle between the B° momentum vector and the vector between the 

92 primary pp interaction and the B° decay vertex is small. Candidates are retained if their 

93 K + n~ invariant mass is in the range 792 < m(K + 7T~) < 992MeV/c 2 . 

94 A multivariate selection, using a boosted decision tree (BDT) [19) with the AdaBoost 

95 algorithm J20], is applied to further reduce the level of combinatorial background. The BDT 

96 is identical to that described in Ref. |6] . It has been trained on a data sample, corresponding 

97 to 36pb _1 of integrated luminosity, collected by the LHCb experiment in 2010. A sample 

98 of B° — )■ K*°J/%1) (J/tp — > candidates is used to represent the B°—> K*°fi + fi~ signal 

99 in the BDT training. Candidates from the B° — > K*°n + fi^ upper mass sideband (5350 < 

100 m(fT + /T~ fi + < 5600MeV/c 2 ) are used as a background sample. Candidates with 

101 invariant masses below the nominal B° mass contain a significant contribution from partially 

102 reconstructed B decays and are not used in the BDT training or in the subsequent analysis. 

103 They are removed by requiring that candidates have m(K + 7r~ fi + > 5150MeV/c 2 . 

104 The BDT uses predominantly geometric variables, including the variables used in the 

105 pre-selection (above). It also includes information on the quality of the B° vertex and the 
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Figure 1: Distribution of /^ + /U~ versus K + 7r~ fi + fi~ invariant mass of selected B°—> K*°fi + fi~ 
candidates. The vertical lines indicate a ±50MeV/c 2 signal mass window around the nominal 
B° mass. The horizontal lines indicate the two veto regions that are used to remove J/ip and 
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• decays. The B° 

windows. 



K fi + \i signal is clearly visible outside of the J/ip and 



we fit x 2 °f the four tracks. Finally the BDT includes information from the RICH and muon 

107 systems on the likelihood that the kaon, pion and muons are correctly identified. Care has 

108 been taken, using simulated events, to ensure that the BDT does not preferentially select 

109 regions of q 2 , K + n~ fi + jjT invariant mass or of the K + n~ fi + jjT angular distribution. The 
no multivariate selection removes 88% of the background that remains after the pre-selection 
in and retains 78% of the signal. 

112 Figure [T] shows the versus K + ix~ fi + fi~ invariant mass of the selected candidates, 

in The B°— > K*°fi + fi~ signal, which peaks in K + ir~fi + fi~ invariant mass, and populates the 

lu full range of the dimuon invariant mass range, is clearly visible. 



us 4 Exclusive and partially reconstructed backgrounds 

lie Several sources of peaking background have been studied using samples of simulated 

117 events, corrected to reflect the difference in particle identification (and misidentification) 

us performance between the data and simulation. Sources of background that are not reduced 

n9 to a negligible level by the pre- and multivariate-selections are described below. 

120 The decays B° -»■ K*°J/tfj and B° -»■ K*°ip(2S), where J/0 and ip(2S) -»■ 

121 are removed by rejecting candidates with 2946 < m(p + fi~) < 3176MeV/c 2 and 
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3586 < m(/i + /i _ ) < 3766MeV/c 2 . These vetoes are extended downwards by 150MeV/c 2 
in m(/i + yU _ ) for B° — ► K*°fi + fi~ candidates with masses 5150 < m(K + Tr~ /i + < 
5230MeV/c 2 to account for the radiative tails of the J/0 and ip(2S) mesons. They are 
also extended upwards by 25MeV/c 2 for candidates with masses above the B° mass to 
account for the small percentage of J/0 or ip(2S) decays that are misreconstructed at 
higher masses. The decay B° — > K* J/ip is used throughout this analysis as a control 
channel. The J/ip and ip(2S) vetoes are shown in Fig. [TJ 

The decay B° — » K*°J/tp can also form a source of peaking background if the kaon or 
pion is misidentified as a muon and swapped with one of the muons from the J/0 decay. 
This background is removed by rejecting candidates that have a K + fi~ or 7i~fi + invariant 
mass (where the kaon or pion is assigned the muon mass) in the range 3036 < m(fi + fi~) < 
3156MeV/c 2 if the kaon or pion can also be matched to hits in the muon stations. A 
similar veto is applied for the decay B° — > K*°ip(2S). 

The decay B®—} <fifi + fi~ , where 0— > K + K~~ , is removed by rejecting candidates if the 
K + 7i~ mass is consistent with originating from a 0— >■ K + K~ decay and the pion is kaon-like 
according to the RICH detectors. A similar veto is applied to remove A® — > A* (1520) fi + fi~ 
(A* (1520) ^ pK~) decays. 

There is also a source of background from the decay B + — > K + fi + fi~ that appears 
in the upper mass sideband and has a peaking structure in cosOk- This background 
arises when a K*° candidate is formed using a pion from the other B decay in the event 
and is removed by vetoing events that have a K + fi + fi~ invariant mass in the range 
5230 < m(K + fi + fi~) < 5330 MeV/c 2 . The fraction of combinatorial background candidates 
removed by this veto is small. 

After these selection requirements the dominant sources of peaking background are 
expected to be from the decays B° — > K*°J/ip (where the kaon or pion is misidentified 
as a muon and a muon as a pion or kaon), B® — > and B® — > K*°fi + fi~ at the 

levels of (0.3 ± 0.1)%, (1.2 ± 0.5)% and (1.0 ± 1.0)%, respectively. The rate of the decay 



B° s — > K*°[i + fi~ is estimated using the fragmentation fraction f s /fd 21 and assuming the 
branching fraction of this decay is suppressed by the ratio of CKM elements |Vtd/Vf S | 2 
with respect to B° — > K*°fi + fi~ . The full size of the effect is assigned as a systematic 
uncertainty. Finally, the probability for a decay B° — > K*°n + fi~ to be misidentified as 
B°—> K*°fi + ii~ is estimated to be (0.85 ± 0.02)% using simulated events. 



154 5 Detector acceptance and selection biases 

155 The geometrical acceptance of the detector, the trigger, the event reconstruction and 

156 selection can all bias the angular distribution of the selected candidates. At low q 2 there 

157 are large distortions of the angular distribution at extreme values of cos#£ (| cos^| ~ 1). 

158 These arise from the requirement that muons have momentum p > 3 GeV/ c to reach the 

159 LHCb muon system. Distortions are also visible in the cos 9k angular distribution. They 
wo arise from the momentum needed for a track to reach the tracking system downstream of 
lei the dipole magnet, and from the impact parameter requirements made in the pre-selection. 
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162 The acceptance in cos 9 k is asymmetric due to the momentum imbalance between the 

163 pion and kaon from the K*° decay in the laboratory frame. 

164 Acceptance effects are accounted for, in a model-independent way, in the angular 

165 and differential branching fraction analyses by weighting candidates by the inverse of 
lee their efficiency determined from simulation. This takes into account the variation of the 

167 acceptance in q 2 to give an unbiased estimate of the observables over the q 2 bin. The 

168 candidate weights are normalised such that they have mean 1.0. The resulting distribution 

169 of weights in each q 2 bin has a root-mean-square in the range 0.2 — 0.4. Only 2% of 

170 candidates have weights larger than 2.0. 

171 The weights are determined using a large sample of simulated three-body B° — > 

172 K*°fi + fi~ phase-space decays. They are determined separately in fine bins of q 2 with 

173 widths: 0.1 GeV 2 /c 4 for q 2 < 1 GeV 2 /c 4 ; 0.2 GeV 2 /c 4 in the range 1< q 2 < 6 GeV 2 /c 4 ; and 

174 0.5GeV 2 /c 4 for q 2 > 6GeV 2 /c 4 . Inside the q 2 bins, the angular acceptance is assumed 

175 to factorise such that s(cos 9^, cos 9 k,4>) = £(cos 9t)s(cos 9k)^{4 > )- This factorisation is 

176 validated at the level of 5% in the phase-space sample. The treatment of the event weights 

177 is discussed in more detail in Sec. 

178 Event weights are also used to account for the fraction of background candidates 

179 that were removed in the lower mass (m(K + ir~ [i + < 5230MeV/c 2 ) and upper mass 
wo (m(K + ir~ fi + > 5330 MeV/c 2 ) sidebands by the J/ip and ip(2S) vetoes described in 

181 Sec. [I] (and shown in Fig. [I]). In each q 2 bin, a linear extrapolation in q 2 is used to estimate 

182 this fraction and the resulting event weights. The weights are largest for the q 2 bin between 

183 the J/iJj and ip(2S) resonances. 



184 6 Differential branching fraction 

185 The angular and differential branching fraction analyses are performed in six bins of q 2 , 

186 which are the same as those used in Ref. |8j. The K + tt~ fi + invariant mass of candidates 

187 in these q 2 bins is shown in Fig. [2] 

188 The number of signal candidates in each of the q 2 bins is estimated by performing an 

189 extended unbinned maximum likelihood fit to the K + tt~ fi + fi~ invariant mass distribution. 

190 The signal shape is taken from a fit to the B° — > K*°J/ip control sample and is parameterised 



200 
201 



191 by the sum of two Crystal Ball 22 functions that differ only by the width of the Gaussian 

192 component. The combinatorial background is described by an exponential distribution. 

193 The decay -6°—)- K*°n + [i~ , which forms a peaking background, is assumed to have a shape 

194 identical to that of the B° — > K*°fi + fi~ signal, but shifted in mass by the B° s — B° mass 

195 difference [23]. Contributions from the decays B® — > (f)[i + ii~ and B° — > K*°J/ip (where 
we the fi~ is swapped with the ir~) are also included. The shapes of these backgrounds are 

197 taken from samples of simulated events. The sizes of the B® — >■ K*°fi + fi~, B®— > <f)fi + fi~ 

198 and B°— > K*°J/ip backgrounds are fixed with respect to the B°— > K*°fi + fi~ signal in the 

199 likelihood fit according to the ratios described in Sec. |3J These backgrounds are varied to 
evaluate the corresponding systematic uncertainty. The resulting signal yields are given in 
Table [l] In the full 0.1 < q 2 < 19.0 GeV 2 /c 4 range, the fit yields 883 ± 34 signal decays. 
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Figure 2: Invariant mass distributions of i^ + 7r~/i + /i~ candidates in the six g 2 bins used in the 
analysis. The candidates have been weighted to account for the detector acceptance (see text). 
Contributions from exclusive (peaking) backgrounds are negligible after applying the vetoes 
described in Sec. HJ 



202 The differential branching fraction of the decay B° — > K*°{i + fi~ , in each q 2 bin, is 

203 estimated by normalising the B° — > K*°fi + fi~ yield, N sig , to the total event yield of the 

204 B° — > K* J/ip control sample, Nk*°m, and correcting for the relative efficiency between 

205 the two decays, £k*°j/iI> / ^ k*° ^+ ^- , 

^ = 2 1 2 / Sig l K ** m x B(B^ K*°J/ip) x B{JM> -> fyr) . (3) 
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Table 1: Signal yield (iV S i g ) and differential branching fraction (dB/dq 2 ) of the B°—> K*°^ + fi" 
decay in the six q 2 bins used in this analysis. Results are also presented in the 1 < q 2 < 6 GeV 2 /c 4 
range where theoretical uncertainties are best controlled. The first and second uncertainties are 
statistical and systematic. The third uncertainty comes from the uncertainty on the B°—> K*°J/tp 
and J/ip — > fi + \i~ branching fractions. The final uncertainty on dB/dq 2 comes from an estimate 
of the pollution from non-K*° B° -> K + -K~n + pT decays in the 792 < m(K + Tr~) < 992MeV/c 2 
mass window (see Sec. |7.3.1). 



q 2 ( GeV 2 /c 4 ) 


N ■ 


dB/dq 2 (10- 7 GeV~ V) 


0.10- 


2.00 


140 ± 13 


0.60 ±0.06 ±0.05 ±0.04 


2.00- 


4.30 


73 ±11 


0.30 ±0.03 ±0.03 ±0.02 


4.30- 


8.68 


271 ± 19 


0.49 ±0.04 ±0.04 ±0.03 


10.09- 


12.86 


168 ± 15 


0.43 ±0.04 ±0.04 ±0.03 t° f 3 


14.18 - 


16.00 


115 ± 12 


0.56 ±0.06 ±0.04 ±0.04 


16.00- 


19.00 


116 ± 13 


0.41 ±0.04 ±0.04 ±0.03 


1.00- 


6.00 


197 ±17 


0.34 ±0.03 ±0.04 ±0.02 
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214 
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The branching fractions B(B° ->■ K*°J/ip) and B(J/ip ->■ are (1.31 ± 0.03 ± 0.08) x 



207 10~ 3 (24J and (5.93 ± 0.06) x 10~ 2 |23j, respectively. 

208 The efficiency ratio, £k*°m /^'V+m - ' depends on the unknown angular distribution of 

209 the B° —)• K*°fi + fi~ decay. To avoid making any assumption on the angular distribution, 

210 the event-by-event weights described in Sec. [5] are used to estimate the average efficiency 
2n of the B° — > K* J/ip candidates and the signal candidates in each q 2 bin. 

212 The resulting differential branching fraction of the decay B° — > K*°fi + fi~ is shown in 

213 Fig. [3] and in Table [1} The bands shown in Fig. [3] indicate the theoretical prediction for 
the differential branching fraction. The calculation of the bands is described in Ref. 25 .^] 



215 In the low q 2 region, the calculations are based on QCD factorisation and soft collinear 

216 effective theory (SCET) |27|, which profit from having a heavy B° meson and an energetic 

217 K*° meson. In the soft-recoil, high q 2 region, an operator product expansion (OPE) in 

218 inverse 6-quark mass (l/m&) and 1/ ^/q 2 is used to estimate the long-distance contributions 

219 from quark loops |28|[29] . No theory prediction is included in the region close to the 
narrow cc resonances (the J/ip and ip(2S)) where the assumptions from QCD factorisation, 
SCET and the OPE break down. The form- factor calculations are taken from Ref. [30~] . 
A dimensional estimate is made of the uncertainty on the decay amplitudes from QCD 

223 factorisation and SCET of (D(AQ CD /mb) [31]. Contributions from light-quark resonances 

224 at large recoil (low q 2 ) have been neglected. A discussion of these contributions can be 

225 found in Ref. (32l. 



2 A consistent set of SM predictions, averaged over each q 2 bin, have recently also been provided by the 
authors of Ref. 26 



8 




Figure 3: Differential branching fraction of the B° — > K*°fi + fj,~ decay as a function of the 
dimuon invariant mass squared. The data are overlaid with a SM prediction (see text) for the 
decay (light-blue band). A rate average of the SM prediction across each q 2 bin is indicated by 
the dark (purple) rectangular regions. No SM prediction is included in the region close to the 
narrow cc resonances. 



6.1 Systematic uncertainties on the differential branching frac- 
tion 

The largest sources of systematic uncertainty on the differential branching fraction come 
from the ~ 6% uncertainty on the combined B° — > K*°J/ip and J/ip — > fi + fi~ branching 
fractions and from the uncertainty on the pollution of non- K*° decays in the 792 < 



m{K + -K^) < 992MeV/c 2 mass window (see Sec. 7.3.1). No S-wave correction needs to 
be applied to the observed yield of B° — >■ K*°J/i/j decays because the branching fraction 
used in the normalisation (from Ref. |24|) corresponds to a measurement of the decay 
B° — > K + 7i~J/if) over the same mass window used in this analysis. 

The uncertainty associated with the corrections used to calibrate the performance 
of the simulation is 1 — 2%. Varying the level of the peaking backgrounds within their 
uncertainties changes the differential branching fraction by 1%. In the simulation a small 
variation in the K + tt~ fi + invariant mass resolution is seen between B°— > K*°J/ifj and 
B° — > K*°fi + fi~ decays at low or high q 2 . Conservatively, the width of the signal mass 
model is varied by up to 5% in every q 2 bin and the deviation from the nominal fit is 
assigned as a systematic uncertainty. 
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242 7 Angular analysis 



243 The B° — >■ angular distribution is reduced to 



d 4 r 



dr/dg 2 d g2 d cos 0£ d cos Ok 1 



9 

16/T 



F L cos 2 ^ + ^(l-F L )(l-cos 2 ^; 



(4) 



F L cos 2 ^ x (2cos 2 ^ - 1) + 

I(l-F L )(l-cos 2 0£)(2co S 2 ^-l) 

S 3 (l -cos 2 9 K )(l - cos 2 ^) cos 20 + 
4 

-y4 FB (l — cos 2 0^-) cos 9i + 
3 

A 9 (l - cos 2 9 K ) (1 - cos 2 £ ) sin 20 

244 after requiring that g 2 ^> 4m 2 and applying the transformation from <fi — y <fi described by 

245 Eq. [2} The parameters A F b, F^, S3 and Ag must satisfy the constraints 



\A FB \ < -(1 - F h ) 



F h ) and \S 3 \ < ^(1 



246 for Eq. [4] to remain positive in all regions of the allowed phase space. These relationships 

247 are automatically taken into account in the transformations 
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A FB = -(1 - F L )A* G and S 3 =-(1-F L )A 

which map A F b and S3 to the theoretically cleaner transverse observables, A^ e and A 2 -,. 
There are less trivial relationships between 5*3, A 9 and A FB , which all depend on the 
same underlying decay amplitudes (see for example Ref. [33] ). These relationships are not 
reproduced here. 

In each of the q 2 bins, A F b (A^ e ), ^l, S3 (A 2 -) and Ag are estimated by performing 
an unbinned maximum likelihood fit to the cos6^, cos6x and distributions of the 
B° — > K*°fi + fi~ candidates. The K + n~ fi + fi~ invariant mass of the candidates is also 
included in the fit to separate between signal- and background-like candidates. The 
background angular distribution is described using the product of three second-order 
Chebychev polynomials under the assumption that the background can be factorised into 
three single angle distributions. This assumption has been validated on the data sidebands 
(5350 < m(K + n-fi + fi-) < 5600MeV/c 2 ). A dilution factor (V = 1 - 2u) is included in 
the likelihood fit for A f -q and Ag, to account at first order for the small probability (lS) for 
a decay B° — y K*°fi + fi~ to be misidentified as B° — y K*°fi + [i~ . The value of u is fixed to 
0.85% in the fit (see Sec.|§. 

Two fits to the dataset are performed: one, with the signal angular distribution 
described by Eq. |4j to measure F L , A FB , S 3 and A 9 and a second replacing A FB and S3 
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265 with the observables A^ e and A\. The angular observables vary with q 2 within the q 2 bins 

266 used in the analysis. The measured quantities therefore correspond to averages over these 

267 q 2 bins. For the transverse observables, where the observable appears alongside 1 — in 

268 the angular distribution, the averaging is complicated by the q 2 dependence of both the 

269 observable and F^. In this case, the measured quantity corresponds to a weighted average 

270 of the transverse observable over q 2 , with a weight (1 — F^dT/dq 2 . 

in 7.1 Statistical uncertainty on the angular observables 

272 The results of the angular fits are presented in Tableland in Figs. H and |5j The 68% 

273 confidence intervals are estimated using pseudo-experiments and the Feldman-Cousins 

274 technique |34|{^] This avoids any potential bias on the parameter uncertainty that could 

275 have otherwise come from using event weights in the likelihood fit or from boundary 

276 issues arising in the fitting. The observables are each treated separately in this procedure. 

277 For example, when determining the interval on Afb, the observables F^, S3 and Ag are 

278 treated as if they were nuisance parameters. At each value of the angular observable 

279 being considered, the maximum likelihood estimate of the nuisance parameters (which also 

280 include the background parameters) is used when generating the pseudo-experiments. The 

281 resulting confidence intervals do not express correlations between the different observables. 

282 The treatment of systematic uncertainties on the angular observables is described in 

283 Sec. Ol 

284 The final column of Table [2] contains the p- value of the SM point in each q 2 bin, which 

285 is defined as the probability to observe a difference between the log-likelihood of the SM 

286 point compared to the best fit point larger than that seen in the data. They are estimated 

287 in a similar way to the Feldman-Cousins intervals by: generating a large ensemble of 

288 pseudo-experiments, with all of the angular observables fixed to the central value of the SM 

289 prediction; and performing two fits to the pseudo-experiments, one with all of the angular 

290 observables fixed to their SM values and one varying them freely. The data are then 

291 fitted in a similar manner and the p- value estimated by comparing the ratio of likelihoods 

292 obtained for the data to those of the pseudo-experiments. The p- values lie in the range 

293 0.18 — 0.72 and indicate good agreement with the SM hypothesis. 

294 As a cross-check, a third fit is also performed in which the sign of the angle (f> for B° 

295 decays is flipped to measure Sg in place of Ag in the angular distribution. The term 5*9 is 

296 expected to be suppressed by the size of the strong phases and be close to zero in every q 2 

297 bin. Aps has also been cross-checked by performing a counting experiment in bins of q 2 . 

298 A consistent result is obtained in every bin. 



3 Nuisance parameters are treated according to the "plug-in" method (see, for example, Ref. 35 ). 
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Table 2: Fraction of longitudinal polarisation of the K*°, F^, dimuon system forward-backward 
asymmetry, Ap-Q and the angular observables S3, Sg and Ag from the B° — > K*°/j, + fi~ decay in 
the six bins of dimuon invariant mass squared, q 2 , used in the analysis. The lower table includes 
the transverse observables A^ e and AL, which have reduced form- factor uncertainties. Results 
are also presented in the 1 < q 2 < 6GeV 2 /c 4 range where theoretical uncertainties are best 
controlled. In the large-recoil bin, 0.1 < q 2 < 2.0GeV 2 /c 4 , two results are given to highlight the 
size of the correction needed to account for changes in the angular distribution that occur when 



n 07+0.10 
u -°' -0.09 



-0.04 
-0.03 



q 2 < lGeV 2 /c 4 (see Sec. 7.2). The final column contains the p- value for the SM point (see text). 

No SM prediction, and consequently no p- value, is available for the 10.09 < q 2 < 12.86 GeV 2 /c 4 
range. 

q 2 ( GeV 2 /c 4 ) /•', .4 
0.10-2.00 
(uncorrected) 
0.10-2.00 
(corrected) 

2.00 -4.30 0.74 ±g;^° ±°;g| -0.20 

4.30 - 8.68 0.57 t^j^ 0.16 

10.09 - 12.86 0.48l£°J 0.28 

14.18 - 16.00 0.33 ±g:g| ±g;g| 0.51 

16.00 - 19.00 0.38 ±8:87 ±8:83 0-30 

1.00 - 6.00 0.65 ±S-2| ±S-2| -0.17 



n 07+0.10 +0.04 

u -°' -0.09 -0.03 



FB 

n 09 +0-12 +0.01 

-U.UZ _Q J2 _Q Q1 



n n9 +0.13 +0.01 
-u.uz _ 0J3 _ _ 01 



1-0.08 
-0.08 
^0.06 
-0.05 
hO.07 
-0.06 
^0.07 
-0.05 
1-0.08 
-0.08 
1-0.06 
-0.06 



+0.01 
-0.01 
+0.01 
-0.01 
+0.02 
-0.02 
+0.02 
-0.02 
+0.01 
-0.02 
+0.01 
-0.01 



^3 

n n A +0.10 +0.01 
-U.U^± _ .io -0.01 



n nc +0.12 +0.01 
-U.UO _ 0J2 -0.01 



-0.04 

0.08 
-0.16 

0.03 
-0.22 

0.03 



+0.10 
-0.06 
+0.07 
-0.06 
+0.11 
-0.07 
+0.09 
-0.10 
+0.10 
-0.09 
+0.07 
-0.07 



+0.01 
-0.01 
+0.01 
-0.01 
+0.01 
-0.01 
+0.01 
-0.01 
+0.02 
-0.01 
+0.01 
-0.01 



Sg 

n n^+ - 10 + - 01 

u - uo -0.09 -0.01 



n« +o i2 +0.01 

u - uu -0.12 -0.01 



-0.03 
0.01 

-0.01 
0.00 
0.06 
0.07 



+0.11 
-0.04 
+0.07 
-0.08 
+0.10 
-0.11 
+0.09 
-0.08 
+0.11 
-0.10 
+0.09 
-0.08 



+0.01 
-0.01 
+0.01 
-0.01 
+0.01 
-0.01 
+0.01 
-0.01 
+0.01 
-0.01 
+0.01 
-0.01 



q 2 ( GeV 2 /c 4 ) 


Ag 




A Re 


p- value 


0.10-2.00 


0.12 


+0.09 +0.01 
-0.09 -0.01 


-0.14 


+0.34 +0.02 
-0.30 -0.02 


-0.04 


+0.26 +0.02 
-0.24 -0.01 


0.18 


(uncorrected) 
















0.10-2.00 


0.14 


+0.11 +0.01 
-0.11 -0.01 


-0.19 


+0.40 +0.02 
-0.35 -0.02 


-0.06 


+0.29 +0.02 
-0.27 -0.01 




(corrected) 
















2.00-4.30 


0.06 


+0.12 +0.01 
-0.08 -0.01 


-0.29 


+0.65 +0.02 
-0.46 -0.03 


-1.00 


+0.13 +0.04 
-0.00 -0.00 


0.57 


4.30 - 8.68 


-0.13 


+0.07 +0.01 
-0.07 -0.01 


0.36 


+0.30 +0.03 
-0.31 -0.03 


0.50 


+0.16 +0.01 
-0.14 -0.03 


0.71 


10.09-12.86 


0.00 


+0.11 +0.01 
-0.11 -0.01 


-0.60 


+0.42 +0.05 
-0.27 -0.02 


0.71 


+0.15 +0.01 
-0.15 -0.03 




14.18 - 16.00 


-0.06 


+0.11 +0.01 
-0.08 -0.01 


0.07 


+0.26 +0.02 
-0.28 -0.02 


1.00 


+0.00 +0.00 
-0.05 -0.02 


0.38 


16.00- 19.00 


0.00 


+0.11 +0.01 
-0.10 -0.01 


-0.71 


+0.35 +0.06 
-0.26 -0.04 


0.64 


+0.15 +0.01 
-0.15 -0.02 


0.28 


1.00-6.00 


0.03 


+0.08 +0.01 
-0.08 -0.01 


0.15 


+0.39 +0.03 
-0.41 -0.03 


-0.66 


+0.24 +0.04 
-0.22 -0.01 


0.72 
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Figure 4: Fraction of longitudinal polarisation of the K*°, Fl, dimuon system forward-backward 
asymmetry, Afb and the angular observables S3 and Ag from the B° — > K*°fj, + fi~ decay as a 
function of the dimuon invariant mass squared, q 2 . The lowest q 2 bin has been corrected for 
the threshold behaviour described in Sec. |7.2| The experimental data points overlay the SM 
prediction described in the text. A rate average of the SM prediction across each q 2 bin is 
indicated by the dark (purple) rectangular regions. No theory prediction is included for Ag, 
which is vanishingly small in the SM. 



299 7.2 Angular distribution at large recoil 

300 Equation p] assumes that q 2 4m 2 , where is the muon mass. This assumption is 

301 valid for (p > lGeV 2 /c 4 but breaks down in the largest recoil, lowest q 2 bin. In the 

302 0.1 < q 2 < 2.0GeV 2 /c 4 range, the angular terms contain an additional q 2 dependence [l], 

303 proportional to 

1 ~ 4m 2 /g 2 (l-4m 2 /g 2 )V 2 

1 + 2m 2 l /q 2 ° r 1 + Imljq 2 1 ' 

304 depending on the angular term. 

305 As q 2 tends to zero, these threshold terms dilute the influence of the observables on 

306 the angular distribution. This dilution leads to a discrepancy between the experimental 

307 definition of the observables, where this additional q 2 dependence is neglected, and the 
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Figure 5: Transverse asymmetries A 2 -, and A^ e as a function of the dimuon invariant mass 
squared, q 2 , in the B° — > K*°/j, + fi~ decay. The lowest q 2 bin has been corrected for the threshold 



behaviour described in Sec. 7.2 The experimental data points overlay the SM prediction that is 
described in the text. A rate average of the SM prediction across each q 2 bin is indicated by the 
dark (purple) rectangular regions. 



308 definition of the observables in terms of the K*° transversity amplitudes. The fraction of 

309 longitudinal polarisation of the K*° meson, F^, is the only observable that is unaffected 

310 by the additional terms; sensitivity to arises mainly through the shape of the cos 9 k 

311 distribution and this shape remains the same whether the threshold terms are included or 

312 not. 

313 In order to estimate the size of the dilution, it is assumed that A$ and A\ are constant 

314 over the 0.1 < q 2 < 2GeV 2 /c 4 region and A^ c rises linearly (with the constraint that 

315 A^ e = at q 2 = 0). The dilution of S3 (Afb) is derived from the relationship between 

316 the observable and A\ {A^f\ Even though F^ is in itself unbiased, an assumption needs 

317 to be made about the q 2 dependence of F L when determining the dilution of the other 

318 observables. An empirical model, 

FL(q 2 ) = , (6) 

1 + aq z 

319 is used. This functional form displays the correct behaviour since it tends to zero as 

320 q 2 — > and rises slowly over the q 2 bin. The coefficient a = 0.67 1° 30 ^ s estimated from 

321 the average of F h obtained from the fit to the 0.1 < q 2 < 2.0 GeV 2 /c 4 bin. 

322 To remain model independent, the dilution factors are estimated by summing over 

323 the (background subtracted) signal candidates in the 0.1 < q 2 < 2GeV 2 /c 4 region. An 

324 example of how this is performed is given in Appendix |B~| for the observable A\. The 

325 typical size of the correction is 10 — 20%. The values of the angular observables, after 

326 correcting for the dilution, are included in Table [2] A similar factor is also applied to 

327 the statistical uncertainty on the fit parameters to scale them accordingly. No systematic 

328 uncertainty is assigned to this correction. 

329 The procedure to calculate the dilution factors has been validated using large samples 
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330 of simulated events, generated according to the SM prediction and several other scenarios 

331 in which large deviations from the SM expectation of the angular observables are possible. 

332 In all cases an unbiased estimate of the angular observables is obtained after applying the 

333 correction procedure. Different hypotheses for the q 2 dependence of Fl, Afb and do 

334 not give large variations in the size of the correction factors. 



7.3 Systematic uncertainties in the angular analysis 

Sources of systematic uncertainty are considered if they introduce either an angular or 
q 2 dependent bias to the acceptance correction. Moreover, three assumptions have been 
made that may affect the interpretation of the result of the fit to the K + n~ fi + fi~ invariant 
mass or angular distribution: that q 2 ^> 4m 2 ; that there are equal numbers of B° and 
B° decays; and that there is no contribution from non-i^*° B° — > K + tt~ decays in 
the 792 < m(K + 7T~) < 992MeV/c 2 mass window. The first assumption was addressed in 



Sec. 1_J2_ and no systematic uncertainty is assigned to this correction. The number of B° 
and B° candidates in the data set is very similar [36] . The resulting systematic uncertainty 



is addressed in Sec. |7.3.2[ The final assumption is discussed in Sec. |7.3.1| below. 

The full fitting procedure has been tested on B°— >■ K*°J/tp decays. In this larger data 
sample, A FB is found to be consistent with zero (as expected) and the other observables 
are in agreement with the results of Ref. (37). There is however a small discrepancy 
between the expected parabolic shape of the cos 9k distribution and the distribution of 
the B° — > K*°J/ip candidates after weighting the candidates to correct for the detector 
acceptance. This percent-level discrepancy could point to a bias in the acceptance model. 
To account for this discrepancy, and any breakdown in the assumption that the efficiencies 
in cos6^, cosOk and <p are independent, systematic variations of the weights are tried in 
which they are conservatively rescaled by 10% at the edges of cos^, cos 6^ and with 
respect to the centre. Several possible variations are explored, including variations that 
are non-factorisable. The variation which has the largest effect on each of the angular 
observables is assigned as a systematic uncertainty. The resulting systematic uncertainties 
are at the level of 0.01 — 0.03 and are largest for the transverse observables. 

The uncertainties on the signal mass model have little effect on the angular observables. 
Of more importance are potential sources of uncertainty on the background shape. In the 
angular fit the background is modelled as the product of three second-order polynomials, 
the parameters of which are allowed to vary freely in the likelihood fit. This model 
describes the data well in the sidebands. As a cross-check, alternative fits are performed 
both using higher order polynomials and by fixing the shape of the background to be flat 
in cosOe, cos 6k and (p. The largest shifts in the angular observables occur for the flat 
background model and are at the level of 0.01 — 0.06 and 0.02 — 0.25 for the transverse 
observables (they are at most 65% of the statistical uncertainty). These variations are 
extreme modifications of the background model and are not considered further as sources 
of systematic uncertainty. 

The angular distributions of the decays B° s — > and B° s — > K*°ji + \i~ are both 

poorly known. The decay B° s — > K*°[i + fi~ is yet to be observed. A first measurement of 
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Table 3: Systematic contributions to the angular observables. The values given are the magnitude 
of the maximum contribution from each source of systematic uncertainty, taken across the six 
principal q 2 bins used in the analysis. 



Source 


^4fb 


F h 


s 3 


<$9 


A 9 




A Re 


Acceptance model 


0.02 


0.03 


0.01 


< 0.01 


< 0.01 


0.02 


0.01 


Mass model 


< 0.01 


< 0.01 


< 0.01 


< 0.01 


< 0.01 


< 0.01 


< 0.01 


B° -> B° mis-id 


< 0.01 


< 0.01 


< 0.01 


< 0.01 


0.01 


< 0.01 


< 0.01 


Data-simulation diff. 


0.01 


0.03 


0.01 


< 0.01 


< 0.01 


0.03 


0.01 


Kinematic reweighting 


< 0.01 


0.01 


< 0.01 


< 0.01 


< 0.01 


0.01 


< 0.01 


Peaking backgrounds 


0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


S-wave 


0.01 


0.01 


0.02 


0.01 


< 0.01 


0.05 


0.04 


B°-B° asymmetries 


< 0.01 


< 0.01 


< 0.01 


< 0.01 


< 0.01 


< 0.01 


< 0.01 



371 



Bg — > cj)fj, + fi~ has been made in Ref. 38 . In the likelihood fit to the angular distribution 

372 these backgrounds are neglected. A conservative systematic uncertainty on the angular 

373 observables is assigned at the level of < 0.01 by assuming that the peaking backgrounds 

374 have an identical shape to the signal, but have an angular distribution in which each of 

375 the observables is either maximal or minimal. 

376 Systematic variations are also considered for the data-derived corrections to the 

377 simulated events. For example, the muon identification efficiency, which is derived from 

378 data using a tag-and-probe approach with J/ip decays, is varied within its uncertainty 

379 in opposite direction for high (p > lOGeV/c) and low (p < lOGeV/c) momentum muons. 

380 Similar variations are applied to the other data-derived corrections, yielding a combined 

381 systematic uncertainty at the level of 0.01 — 0.02 on the angular observables. The correction 

382 needed to account for differences between data and simulation in the B° momentum 

383 spectrum is small. If this correction is neglected, the angular observables vary by at most 

384 0.01. This variation is associated as a systematic uncertainty. 

385 The systematic uncertainties arising from the variations of the angular acceptance are 

386 assessed using pseudo-experiments that are generated with one acceptance model and 

387 fitted according to a different model. Consistent results are achieved by varying the event 

388 weights applied to the data and repeating the likelihood fit. 

389 A summary of the different contributions to the total systematic uncertainty can be 

390 found in Table |3j The systematic uncertainty on the angular observables in Table [2] is the 

391 result of adding these contributions in quadrature. 



392 7.3.1 Influence of S-wave interference on the angular distribution 

393 The presence of a non-i^* B° — > K + n~ fi + {i~ component, where the K +/ k~ system is in 

394 an S-wave configuration, modifies Eq. [4] to 
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400 



424 
425 



d 4 r 

= 1-^8 



dF/dg 2 d g 2 d cos 9e d cos Ok d(b 

9 

+ 



i d 4 r 



16n 



Ldr/dg 2 d g 2 d cos Q( d cos Ok 
2 _ o 4 

"scos^k 



F s (l - cos 2 9t) + ^-A s cos Ml - cos 2 6 t 



(7) 



395 where F s is the fraction of B° -> K+tt' [i + n~ S-wave in the 792 < m(F+7r-) < 992 MeV/c 2 

396 window. The partial width, Y' , is the sum of the partial widths for the B° — > K*°fi + [i~ 

397 decay and the B° — > K + ir~n + [i~ S-wave. A forward-backward asymmetry in cos6k, A$, 

398 arises due to the interference between the longitudinal amplitude of the K*° and the 



399 S-wave amplitude 39 



42 



The S-wave is neglected in the results given in Table [2] To estimate the size of the 

401 S-wave component, and the impact it might have on the B°— > K*°fi + [i~ angular analysis, 

402 the phase shift of the K*° Breit-Wigner function around the K*° pole mass is exploited. 

403 Instead of measuring Fg directly, the average value of A$ is measured in two bins of K + n~ 

404 invariant mass, one below and one above the K*° pole mass. If the magnitude and phase 

405 of the S-wave amplitude are assumed to be independent of the K + tt~ invariant mass in 

406 the range 792 < m(K + 7r~) < 992MeV/c 2 , and the P-wave amplitude is modelled by a 

407 Breit-Wigner function, the two A$ values can then be used to determine the real and 

408 imaginary components of the S-wave amplitude (and F s )|^] 

409 For a small S-wave amplitude, the pure S-wave contribution, Fg, to Eq. [7] has only a 

410 small effect on the angular distribution. The magnitude of As arising from the interference 

411 between the S- and P-wave can however still be sizable and this information is exploited 

412 by this phase-shift method. The method, described above, is statistically more precise 

413 than fitting Eq. [7] directly for A$ and F$ as uncorrelated variables. For the B° — > K*°J/i/j 

414 control mode, the gain in statistical precision is approximately a factor of three. 

415 Due to the limited number of signal candidates that are available in each of the 

416 q 2 bins, the bins are merged in order to estimate the S-wave fraction. In the range 

417 0.1 < q 2 < 19 GeV 2 /c 4 , F s = 0.03 ± 0.03, which corresponds to an upper limit of F s < 0.04 

418 at 68% confidence level (CL). The procedure has also been performed in the region 

419 1 < q 2 < 6GeV 2 /c 4 , where both F L and F s are expected to be enhanced. This gives 

420 F s = 0.04 ± 0.04 and an upper limit of F s < 0.07 at 68% CL. Conservatively, F s = 0.07 

421 is used to estimate a systematic uncertainty on the differential branching fraction and 

422 angular analyses. The B°— > K*°J/i/j data has been used to validate the method. 

423 For the differential branching fraction analysis, Fs scales the observed branching 
fraction by up to 7%. For the angular analysis, Fs dilutes A FB , S 3 and A 9 . The impact 
on Fl however, is less easy to disentangle. To assess the possible size of a systematic bias, 



4 In the decay B° —> K*°fi + /i~ there are actually two pairs of amplitudes involved, left- and right-handed 
longitudinal amplitudes and left- and right-handed S-wave amplitudes (where the handedness refers to 
the chirality of the dimuon system). In order to exploit the interference and determine Fg it is assumed 
that the phase difference between the two left-handed amplitudes is the same as the difference between 
the two right-handed amplitudes, as expected from the expression for the amplitudes in Refs. 39 40 . 
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426 pseudo-experiments have been carried out generating with, and fitting without, the S-wave 

427 contribution in the likelihood fit. The typical bias on the angular observables due to the 

428 S-wave is 0.01 — 0.03. 



429 7.3.2 Production, detection and direct CP asymmetries 

430 If the number of B° and B° decays are not equal in the likelihood fit then the terms in 

431 the angular distribution no longer correspond to pure CP averages or asymmetries. They 

432 instead correspond to admixtures of the two, e.g. 

S° bs ^S 3 -A 3 (Acp + kA p + Ad) , (8) 

433 where Aqp is the direct CP asymmetry between B° — > K*°n + and B° — > K*°{i + {i~ 

434 decays, Ap is the production asymmetry between B° and B° mesons, which is diluted by 

435 a factor k due to B° — B° mixing, and An is the detection asymmetry between the B° 

436 and B° decays (which might be non-zero due to differences in the interaction cross-section 

437 with matter between K + and K~ mesons). In practice, the production and detection 

438 asymmetries are small in LHCb and ^4cp is measured to be Ac? = — 0.072±0.040±0.05 [36] , 

439 which is consistent with zero. This, combined with the expected small size of the CP 

440 asymmetry or CP-averaged counterparts of the angular observables measured in this 

441 analysis, reduces any systematic bias to < 0.01 (this is neglected in Table [3]). 



8 Forward- backward asymmetry zero-crossing point 

In the SM, A FB changes sign at a well defined value of q 2 , q^, whose prediction is largely 
free from form-factor uncertainties. It is non-trivial to estimate q 2 from the angular fits 
to the data in the different q 2 bins. Instead, q 2 is estimated by performing a fit to the 
q 2 and K + ir~ invariant mass distribution of the forward-going (cos 6^ > 0) and 
backward-going (cos6^ < 0) candidates. If -Pf(<? 2 ) describes the q 2 dependence of the 
forward-going, and P B {q 2 ) the backward-going candidates then 

, N F P F (q 2 ) - N B P B (q 2 ) 
] ~ N F P F (q 2 ) + N B P B (q 2 ) ■ (9) 

Here, N F and N B are the total number of forward- and backward-going signal decays, 
respectively. The zero-crossing point of A FB is found by solving for the value of q 2 at 
which A FB is zero. 

Using candidates in the range 1.0 < q 2 < 7.8GeV 2 /c 4 and third-order polynomials to 
describe the q 2 dependence of the signal and the background, the zero-crossing point is 
found to be 



q 2 = 4.9 ± 0.9 GeV 2 /c 4 . 

The uncertainty on q 2 , is determined using a bootstrapping technique [43]. The zero- 
crossing point is largely independent of the polynomial order and the q 2 range that is 
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used. This value is cons istent with SM predictions 
3.9-4.4GeV 2 /c 4 " 



4.36 toil GeV 



44 



46 



which are typically in the range 
and have relative uncertainties below the 10% level, example, 
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The systematic uncertainty on the zero-crossing point of the forward-backward asym- 
metry is negligible compared to the statistical uncertainty. To generate a large systematic 
bias, it would be necessary to create an asymmetric acceptance effect in cos 9g that is not 
cancelled when combining B° and B° decays. The combined systematic uncertainty is at 
the level of ±0.05 GeV 2 /c 4 . 



46 5 9 Conclusions 

466 In summary, using a data sample, corresponding to l.Ofb -1 of integrated luminosity, 

467 collected by the LHCb experiment in 2011, the differential branching fraction, dB/dq 2 , of 

468 the decay B° — > K*°fi + fi~ has been measured in bins of q 2 . Measurements of the angular 

469 observables, Afb (^T e )> Fl, 5*3 (A^) and A 9 have also been performed in the same q 2 

470 bins. These are the most precise measurements of dB/dq 2 and the angular observables to 

471 date. A first measurement of the zero-crossing point of the forward-backward asymmetry 

472 has also been performed by exploiting the q 2 distribution of forward- and backward-going 

473 candidates, with the value q 2 . = 4.9 ± 0.9GeV 2 /c 4 . All of the observables are consistent 

474 with SM expectations and together put stringent constraints on the contributions from 

475 new particles to b — > s flavour changing neutral current processes. 
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491 Appendix 



492 A Angular basis 

493 The angular basis used in this paper is illustrated in Fig. |6j The angle 9# is defined as the 

494 angle between the direction of the fi + (fi~) in the dimuon rest frame and the direction 

495 of the dimuon in the B° (B°) rest frame. The angle 8k is defined as the angle between 

496 the direction of the kaon in the K*° (K*°) rest frame and the direction of the K*° (K*°) 

497 in the B° (B°) rest frame. The angle is the angle between the plane containing the fi + 

498 and [A and the plane containing the kaon and pion from the K*°. Explicitly, cos^ and 

499 cos 6^ are defined as 

cos0*-^+ -p^+^J - [p^ ■ [~P B o )) , (10) 

(K*°) ~(B°)\ _ ( ~(K*°) ( ~(K*°) 



cos6 K = >.p^>) = (p£ > ■ [-pf >)) (11) 
for the B° and 

co S «, = • p<f >_) = ■ {-P'^)) , (12) 

oosfe = (fip.f£>) = {pT"> ■ {-V l P)) (13) 
for the B° decay. The definition of the angle is given by 

cos^rfer»xpH.te o )xpf_°))) , (M) 



smi 

for the B° and 



v AB°)\ v (ab°) v A B °)\\ A B °) 



x {P K K + , xP K :- > -P K kJ (15) 



cos0 = ((p^ 0) x pffi) ■ (p^A x ^+ 0) )) , (16) 

sin0 = - ((^ B0) x ^+ 0) J x (pjf_ 0) x p { f+ ] jj ■ pgj (17) 

503 for the 5° decay. The are unit vectors describing the direction of a particle X in 

504 the rest frame of the system Y . In every case the particle momenta are first boosted to 

505 the B° (or B°) rest frame. In this basis, the angular definition for the B° decay is a CP 

506 transformation of that for the B° decay. 
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(a) 9k and 9g definitions for the B° decay 




Figure 6: Graphical representation of the angular basis used for B — > K*°[i + ij,~ and 
B — > K*°fj, + fi~ decays in this paper. The notation h a \> is used to represent the normal to 
the plane containing particles a and b in the B° (or B°) rest frame. An explicit description of 
the angular basis is given in the text. 
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B Angular distribution at large recoil 



508 
509 
510 



An explicit example of the dilution of the angular observables that comes from the threshold 
terms is provided below for A\. Sensitivity to A\ comes through the term in Eq. [I] with 
sin 2 9i sin 2 9^ cos 2(f) angular dependence. In the limit q 2 ^> m 2 this term is simply 



1 



(1 - F L (q 2 )) A%(q 2 ) sin 2 6 t sin 2 9 K cos 2(f) 



5ii and the differential decay width is 



dr 

ay 



\A 



0,L| 



+ L4n. L r + l^j 



+ |A , R | 2 + |A|| iR | 2 + \A ±tK \ 



(18) 



(19) 
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513 
514 



where A , An and Aj_ are the K*° spin-amplitudes and the L/R index refers to the chirality 
of the lepton current (see for example Ref. [lj). If q 2 < 1 GeV 2 /c 4 , these expressions are 
modified to 



1 — 4m 2 /q 2 
2 [1 + 2m 2 /q 2 _ 



(l - F L (q 2 )) A 2 T (q 2 ) sin 2 9 e sin 2 9 K cos 20 



(20) 



and 



dr 

dq 2 



[1 + 2ml/q 2 } (|,4o, L | 2 + |A| |iL | 2 + \A ± , L \ 2 + |A , R | 2 + |A|,, R | 2 + \A ±yR \ 2 ) . (21) 



516 In an infinitesimal window of g 2 , the difference between an experimental measurement 

517 of A\, AL^, in which the threshold terms are neglected and the value of A\ defined in 



5i8 literature is 



a 2 exp 



- 4m 2 jq 2 
_ 1 + 2m 2 /q 2 



(22) 



519 Unfortunately, in a wider q 2 window, the q 2 dependence of F^, A\ and the threshold terms 

520 needs to be considered and it becomes less straightforward to estimate the dilution due to 

521 the threshold terms. If A\ is constant over the q 2 window, 



/•?Lx dr 


"1 — 4m 2 jq 2 ' 


A 2exp J ^2 


1 + 2ml/ q 2 _ 


A\ ~ ril 


lax rl 
li 





[1 - F L (q 2 )} dq 2 



(23) 



522 In practice the integration in Eq. 23 can be replaced by a sum over the signal events in 

523 the q 2 window 
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A 2 exp 



N 

E 

t=0 



l~4m|/g| 
l+2m2/g2 



JV 



(24) 



£(i-*L(fl?))<* 



where is a weight applied to the 2 th candidate to account for the detector and selection 
acceptance and the background in the q 2 window. 

Correction factors for the other observables can be similarly defined if it is assumed that 
they are constant over the q 2 window. In the case of A FB (and A^ e ) that are expected to 
exhibit a strong q 2 dependence, the q 2 dependence of the observable needs to be considered. 
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